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ABSTRACT 


OBJECTIVE: To assess the frequency and patterns of mutations and distribution of alleles in 
Plasmodium. Vivax Dihydrofolate Reductase (DHFR) & Dihydro Pteroate Synthase (DHPS) genes in 
Karachi. 

METHODOLOGY: The descriptive cross-sectional study was carried out at the Department of 
Biotechnology, University of Karachi from June 2014 to February 2016. The sampling technique used 
was Non-probability convenient sampling method. 

A total of n=200 malaria cases irrespective of age and gender from Jinnah Medical College Hospital 
Korangi Karachi were selected. Around 200 yL of blood was collected in micro centrifuge tubes. Among 
them, 36 cases of P. vivax were detected by microscopy. DNA was extracted from these 36 samples 
followed by polymerase chain reaction (PCR) and sequencing of amplified genes of Plasmodium vivax 
Dihydrofolate reductase (PVDHFR) & Plasmodium vivax Dihydro pteroate synthase (PVDHP). 

RESULTS: The analysis of 36 amplified P. vivax isolates revealed presence of various alleles of pvdhfr 
and pvdhps genes producing both wild type (W) and mutant amino acids having single, double & 
multiple mutations at multiple positions. In pydhfr gene, common mutations were observed in ‘S’ amino 
acids at codons 58 &117 and ‘D’ at codons 105 &157. Regarding pvdhps gene, common mutations were 
observed in ‘V’ amino acids at codon 398 and ‘K’ at codon 590. 

CONCLUSION: The frequency and variety of mutations we detected in pvdhfr and pvdhps genes signify 
the drug resistance in plasmodium vivax so there is need of some measures to control drug resistance 
by better formulation using combination of drugs. 
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their own folic acid catalyzed by DHPS and DHFR 
enzymes. Pyrimethamine inhibits the dhfr enzyme and 


INTRODUCTION 


Malaria is highly endemic in developing countries like 
Pakistan’. According to the World malaria report 2018, 
malaria accounted for an estimated 219 million cases 
worldwide as reported in 2017, and there were an 
estimated 435,000 deaths from malaria globall 
Plasmodium vivax is the most widely distributed 
malarial species. Control of P. vivax infection is more 
challenging because of its unique biological and 
epidemiological characteristics. The level of 
parasitemia is low but gametocytes emerge before 
illness manifests and before control measures are 
taken. Also, hypnozoite (dormant liver stage) infection 
occurs out of season and so can extend the 
geographical area affected’. 

P. vivax infection has high morbidity and some 
mortality’. This can be attributed to many factors, 
mainly treatment failure caused by emergence of drug 
resistant strains including antifolate drugs such as 
Sulphadoxine-Pyrimethamine (SP). Persistent use of 
SP produces resistant strains, which develo oP 
sequential mutations in dhfr and dhps genes” 
Unicellular eukaryotes and prokaryotes synthesize 
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sulfadoxine targets the dhps enzyme in the folate 
biosynthetic pathway of the parasite’®. SP is a cheap 
anti malarial drug commonly prescribed against P. 
vivax in Asian countries’. SP resistance in P. vivax is 
reported in many regions of Pakistan'®’*. There are 
several factors which are responsible for the 
emergence of resistance like microscopic misdiagno- 
sis and irrational use of anti malarial drugs especially 
Sulphadoxine, which is also used for the treatment of 
GIT & respiratory tract infections as well as other 
causes of fever and chill'*. Molecular markers for SP 
resistance, mutations in dhfr and dhps genes are 
valuable tools for evaluating SP resistance patterns in 
P. vivax populations’. 

The aim of this study was to evaluate the extent of 
antifolate-resistance in P.vivax species by detecting 
mutations in DHFR and DHPS and to outline the 
changing genetic diversity. 


METHODOLOGY 


The descriptive cross-sectional study was carried out 
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at the Department of Biotechnology, University of 
Karachi. The duration of project carried out was from 
June 2014 to February 2016. The sampling technique 
used was  Non-probability convenient sampling 
method. 

A total n= 200 malaria cases irrespective of age and 
gender were selected from various hospitals & clinics 
of Karachi. Around 200 uL of blood was collected in 
micro centrifuge tubes. Among them, 36 cases of P. 
vivax were detected by microscopic examination of 
Giemsa stained thick and thin smears as described by 
WHO guidelines. 

DNA was extracted from blood samples of these P. 
vivax cases using the QIAGEN DNA extraction kit 
(Germany) According to manufacturers protocol and 
instructions. 

The PCR was performed by using nested PCR for 
confirmation and identification of new genes. The 
confirmation of P. vivax was done by targeting the 18S 
rRNA, as reported by Iqbal et al in 2013"°. Further 
pvdhf and dhps genes were amplified and sequencing 
was performed for detecting mutation or 
polymorphism in these genes. The extracted samples 
of all isolates were then stored at -20°C for further 
use. Quality control of each lot of the QLAamp DNA 
Mini kit is already tested against predetermined 
specification to ensure consistent product quality. 


Amplification of Pvdhfr & Pvdhps Gene 


A total of 36 P. vivax positive cases were selected and 
nested PCR was performed on dhfr and dhps genes. 
The primer sequences are shown in Table I. 


TABLE I: PRIMER SEQUENCES USED FOR 
AMPLIFICATION OF PVDHFR AND PVDHPS 
GENES 


Gene Primer No Primer Sequence 

Pvdhfr 107 F: 5-CACCGCACCAGTTGATTCCT-3’ 
(P) 108 R: 5'-CCTCGGCGTTGTTCTTCT-3’ 
Pvdhfr 109 F:5’-CCCCACCACATAACG AAG-3’ 
(S) 110 R:5’-CCCCACCTTGCTGTAAACC-3’ 
Pvdhps 111 F: 5'-GATGGCGGTTTATTTGTCG-3’ 
(P) 112 R: 5'-GCTGATCTTTGTCTTGACG-3’ 
Pvdhps 113 F: 5'-GCTGTGGAGAGGATGTTC-3’ 
(S) 114 R: 5'-CCGCTCATCAGTCTGCAC-3’ 


The cycling conditions were used _ with slight 
modifications as previously reported by Huang et al in 
2015'°. The PCR was carried out in a tube containing 
50 ul of a reaction mixture made up of the following 
components: 15 pmol of each forward and reverse 
primers for pvdhfr & pvdhps genes as shown in Table 
I. The 500 uM of four deoxynucleotides, 2 U of Taq 
polymerase (Promega), 10 x PCR buffer and 1.5 mM 
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MgClz. The thermal cycler (Master Gradient PCR 
System, Eppendorf AG, Germany) was programmed 
to first incubate the sample for 10 minute at 94°C 
followed by 40 cycles at 94°C for 1 minute, at 45°C for 
2 minutes and at 72°C for 1 minute with final 
extension for 7 minute at 72°C.PCR amplification of 
pvdhfr gene was done at codons 50, 52, 58, 93, 105, 
106, 117,135, 157 and 185 while for pvdhps gene was 
done at positions 389, 391, 397, 398, 399, 400, 401, 
426, 459, 528, 549, 577, 582, 586, 588, 589 and 590. 

The products were identified by electrophoresis on a 
2% agarose gel stained with ethidium bromide, and 
evaluated under transilluminator as described by Ding 
S et al'’. The sizes of PCR amplified product were 
estimated according to the migration pattern of a 1 kb 
DNA ladder (Gibco BRL Life Technologies). The 
amplified products were visualized and photographed 
using DOC gel documentation system (Vilber 
System). The amplified fragment of pvdhfr gene was 
755 bp size shown in Figure | and of pvdhps gene 
was 731 bp shown in Figure Ill. The samples were 
then sent for sequencing. Mutation analysis was done 
by aligning the sequence with the reference 
sequence. 

PVDHFR and PVDHPS Sequence Analysis 

The nested-PCR products of pvdhfr and pvdhps 
genes were subjected to sequence analysis. The 
forward primer was_ purified and sequenced 
commercially by Macrogen Incorporation in Seoul 
Korea. DNA sequences were translated into amino 
acid sequences and compared with the published 
sequences of pvdhfr and pvdhps genes'®. Multiple 
sequence alignment of P.vivax dhfrgene is shown in 
Figure Il, and multiple sequence alignment (MSA) of 
P.vivax dhps gene is shown in Figure VI, V and VI. 


FIGURE I: PYDHFRGEL ELECTROPHORESIS: 
showing fragment size. 

Lanes 1 and 19: band on sample 7 not seen.100 bp 
DNA molecular weight marker 

Lanes 2 — 20 and 26 — band on sample 21 and 26 not 
seen: Amplified pvdhfrgenomic DNA 
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FIGURE Ill: PYVDPS GELECTOPHORESIS: Showing fragment size. 
FIGURE IV: MULTIPLE SEQUENCE ALIGNMENT OF PVDHPS GENE 


Lanes 1 and 25: 100bp DNA molecular weight markers 


FIGURE II: MULTIPLE S 
Lanes 2 


FIGURE V: MULT 
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FIGURE VI: MULTIPLE SEQUENCE ALIGNMENT OF PVDHPS GENE 


‘ 
FAH iA TURAL LAAAEA EAA LTTNG 
: SERUEREAL BASES EL SULLA OOEE 


cAtti 


Bioinformatics tools to observe the pattern of 
resistance in Pvdhfr and Pvdhps Genes 

The DNA sequences were aligned by using various 
web tools. Six frame translation of nucleotide 
sequence into amino acid sequence was done using 
EMBOSS transeq tool of EBI. The suitable frame was 
selected with the help of Blast suite of NCBI. Clustal X 
2.1 was used to perform multiple sequence alignment 
and phylogenetic analysis. Secondary structure 
prediction was done by PSIPRED method, and a 3D 
model was made using Swiss model. Energy calcula- 
tion and Ramachandran plot analysis was done by 
PROSA web server and RAMPAGE, respectively. 


RESULTS 


Of the total 200 cases, only 36 were processed for 
amplification and sequencing. The samples showed 
wild type (W) as well as mutated haplotype of DHFR 
genes when at least one amino acid substitution in the 
related residue was seen. The mutation patterns and 
wild type haplotypes of DHFR gene is summarized in 
Table Il. 


TABLE I: 
MUTATION ANALYSIS OF PVDHFR GENE 
MUTANT TYPE 
SOPON: WILEHIYPE Mutant Type No. of Mutants 
50 N I 2 
52 Vv A 1 
58 S) R 17 
93 S) H 1 
N 18 
105 D Y 1 
H 4 
E 1 
106 K 
Q 1 
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117 S N 20 
135 Vv G 1 
157 D N 21 
185 R K 3 


The most common W type amino acids detected in 


pvdhfr gene were S at codons 58, 93 & 117, D at 
codons 105 & 157 and V at codons 52 & 135. 
Common mutations detected in pvdhfr gene were 
S58R (17/36, 47% of cases), S117N (20/36, 56% of 
cases), D105N (18/36, 50% of cases) and D157N 
(21/36, 58 % of cases). 

In dhfr gene,2 types of mutations were seen at 
position 106 (K106E in 1/36, 2.8 % of cases) (K106Q 
in 1/36, 2.8 % of cases) and 3 types of mutations at 
position 105 (D105N in 18/36, 50 of cases %), (D105Y 
in 1/36, 2.8 % of cases) and (D105H in 4/36, 11.1% of 
cases). 

The mutation patterns and wild type haplotypes 
observed in pvdhps gene are mentioned in Table Ill. 
TABLE Ill: 

MUTATION ANALYSIS OF PVDHPS GENE 


MUTANT TYPE 

CODON WILDTYPE “Witant Type No. of Mutants 
389 N i : 
Ss 1 
301 S T ’ 
397 L Vv i 
398 Vv D eg 
399 M R d 
400 P S d 
401 Vv A dl 
426 A S a 
459 D A i 
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528 
549 
577 
582 
586 
588 
589 
590 K 
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Regarding pvdhps gene, the most common wild type 
amino acids were Gat codons 549, 577 & 589, V at 
codons 398 & 401 and M at codons 399 & 588. In 
pvdhps gene, common mutations were V398D (18/36, 
50% of cases) and K590E (15/36, 42% of cases). 
Codon 389 had 2 types of mutations (N389R in 3/36, 
8.3% of cases) (N389S in 1/36, 2.8% of cases). 


DISCUSSION 


The recent trend of resistance of P. vivax isolates to 
various antimalarials including antifolate drugs like 
Sulphadoxine-Pyrimethamine (SP) is alarming and 
need to be monitored carefully to avoid losing this 
important drug option. Molecular drug resistance 
markers for SP resistance include genes encoding 
DHFR and DHPS enzymes which can be valuable 
tools to serve this purpose as suggested by many 
studies". 

In this study, we examined the distribution and pattern 
of mutant alleles and polymorphisms in pvdhfr and 
pvdhps genes of P. vivax by sequence analysis. The 
established data on P. vivax resistance profile is 
scarce in Pakistan especially in Karachi. Our study 
revealed both distinct polymorphisms of wild type 
allele and also high prevalence of mutant alleles with 
double triple as well as novel mutations in both pvdhfr 
and pfdhps genes. These mutations have shown to 
confer antifolate drug resistance °°. 

We found the most common wild type amino acids in 
pvdhfr gene are S at codons 58, 93 &117, D at 
codons 105 &157 and V at codons 52 & 135. Zakeri S 
20097 in his study observed prevalence of both wild 
type alleles as well as mutants in pvdhfr genes. Auliff 
A et al'®. examined genetic mutations in dhfr and dhps 
in P. vivax samples from six Asian-Pacific countries. 
They observed wild-type dhfr alleles to be present ata 
relatively high proportion in P. vivax parasites from 
China, East Timor, the Philippines, and Vietnam while 
less so from Papua New Guinea and Vanuatu. 

We found high level of mutations in Pvdhfr gene. 
Among mutant strains we found amino acid S and D 
had multiple mutations. Common mutations found in 
Pvdhfr genes were S58R (17/36, 47% cases), S117N 
(20/36, 56% cases), D105N (18/36, 50% cases) and 
D157N (21/36, 58 % cases). Amino acid D showed 3 
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types of mutations at codon105 (D105N, D105Y & 
D105H). 

Many other studies also reported mutations S58R and 
$117NinPvdhfr gene similar to our observations. 
Gresty et al’? reported high prevalence of mutant 
alleles (58R/117N) in pvdhfr genes. Also studies by 
Raza A 2013", Kuesap J 20117’ and Khatoon L 
2009" reported similar mutations in pvdhfr gene. 
Recent studies of Nyunt et al.reported variant of wild 
type and mutant allele withdouble mutations (S58R 
and S117T/N) or quadruple mutations (F57L/I, S58R, 
T61M, and S117T) of pvdhfr gene in Thailand, 
Myanmar, Korea, Cambodia, and India’*. Studies by 
Huang B et al'® found mutations at codons F57I/L, 
S58R, T61M, H99R/S, and S117N/T. 

In contrast to our findings, Khattak AA 20137 
observed low prevalence of SP resistance associated 
mutations in P.vivax in Pakistan. However, there were 
pvdhfr S117N single and S58R plusS117N double 
mutants with highest prevalence in Sindh and Khyber 
Pakhtunkhwa provinces. 

Also studies of Sharifi-Sarasiabi K 2016* from 
Southern Iran and Lekweiry KM et al** from Mauritania 
reported low prevalence of pvdhfr gene mutation 
reflecting reduced antifolate drug resistance and these 
results are contrary to our study. 

Regarding pvdhps gene, our study showed common 
wild type amino acids were Gat codons 549, 577 & 
589, V at codons 398 & 401 and M at codons 399 & 
588. In our study, we observed that most common 
mutations in pvdhps gene were in amino acids V and 
K. Common mutations were V398D (18/36, 50% of 
cases) and K590E (15/36, 42% of cases). A similar 
pattern of significantly higher mutant rate of pvdhps, 
was observed by Nyunt MH et al’ with triple or 
quadruple mutant alleles which accounted for more 
than 50% of all samples. 

Kuesap J 2011*' described mutations at codons 383 
and 553 of Pvdhps gene. Their results suggest that all 
P. vivax isolates in Thailand carry some combination 
of mutations in Pvdhfr and Pvdhps genes. Similar 
results were seen by Pornthanakasem W et al® who 
confirmed the role of key mutations in dhps gene in 
sulfa drug resistance in P. vivax. Enzyme inhibition 
studies with sulfa drugs revealed that differences in 
sulfa drug sensitivity between the wild-type and 
mutant enzymes are responsible for the sulfa drug 
resistance observed for P. vivax. 

In contrast to our study, Gretsy KJ et al'?and Zakeri S 
2009”° observed that all of P.vivax isolates possessed 
a wild type dhps gene and no mutation was seen. 
Similarly, Khattak AA 2013 and Lekweiry KM et al** 
reported that majority of samples contained the 
wild-type haplotype in pvdhps gene. 


CONCLUSION 


The obtained results after multiple sequence 
alignment of pvdhfr gene and changing of amino acid 
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sequence on various positions revealed a large 
variety of mutations in pvdhfr and pvdhps genes. This 
signifies that urgent measures needed for control of 


drug 


resistance and development of effective, 


affordable new anti malarial drugs that should be used 
judiciously. The future of effectual malaria control will 
depend heavily upon appropriate data collection and 
more precise evidence-based decision on changing 
pattern of mutant strain of P. vivax as there is no anti 
malarial vaccine available so far. 
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